Fermionic field theory for trees and forests 
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We prove a generalization of Kirchhoff's matrix-tree theorem in which a large class of combina- 
torial objects are represented by non-Gaussian Grassmann integrals. As a special case, we show 
that unrooted spanning forests, which arise as a q — > limit of the Potts model, can be represented 
by a Grassmann theory involving a Gaussian term and a particular bilocal four-fermion term. We 
show that this latter model can be mapped, to all orders in perturbation theory, onto the A-vector 
model at N = —1 or, equivalently, onto the cr-model taking values in the unit supersphere in R 1 ' 2 . 
It follows that, in two dimensions, this fermionic model is perturbatively asymptotically free. 
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Kirchhoff's matrix-tree theorem 1] and its generaliza- 
tions f|, which express the generating polynomials of 
spanning trees and rooted spanning forests in a graph 
as determinants associated to the graph's Laplacian ma- 
trix, play a central role in electrical circuit theory Q and 
in certain exactly-soluble models in statistical mechanics 
PHi}. Like all determinants, those arising in Kirchhoff's 
theorem can of course be rewritten as Gaussian integrals 
over fermionic (Grassmann) variables. 

In this Letter we prove a generalization of Kirchhoff's 
theorem in which a large class of combinatorial objects 
are represented by suitable non-Gaussian Grassmann in- 
tegrals. Although these integrals can no longer be cal- 
culated in closed form, our identities allow the use of 
field-theoretic methods to shed new light on the critical 
behavior of the underlying geometrical models. 

As a special case, we show that unrooted spanning 
forests, which arise as a q — > limit of the g-state 
Potts model 0, can be represented by a Grassmann 
theory involving a Gaussian term and a particular bilo- 
cal four-fermion term. Furthermore, this latter model 
can be mapped, to all orders in perturbation theory, 
onto the A^-vector model [0(A^)-invariant cr-model] at 
N = — 1 or, equivalently, onto the cr-model taking values 
in the unit supersphere in R 1 ' 2 [OS'P(l|2)-invariant cr- 
model]. It follows that, in two dimensions, this fermionic 
model is perturbatively asymptotically free, in close anal- 
ogy to (large classes of) two-dimensional cr-models and 
four-dimensional nonabelian gauge theories. Indeed, this 
fermionic model may, because of its great simplicity, be 
the most viable candidate for a rigorous nonperturbative 
proof of asymptotic freedom — a goal that has heretofore 
remained elusive in both cr-models and gauge theories. 

The plan of this Letter is as follows: First we prove 
some combinatorial identities involving Grassmann inte- 
grals, culminating in our general formula l|12|) . and show 
how a special case yields unrooted spanning forests. Next 



we show that this latter model can be mapped onto the 
TV- vector model at N = —1, and use this fact to de- 
duce its renormalization-group (RG) flow at weak cou- 
pling. Finally, we conjecture the nonperturbative phase 
diagram in this model. 

Combinatorial Identities. Let G — (V, E) be a finite 
undirected graph with vertex set V and edge set E. As- 
sociate to each edge e a weight w e , which can be a real 
or complex number or, more generally, a formal algebraic 



variable. For i ^ j, let 



be the sum of w e over 



all edges e that connect itoj. The (weighted) Laplacian 
matrix L for the graph G is then defined by = — my 
for i 7^ j, and La = J2k^a w ik- This is a symmetric 
matrix with all row and column sums equal to zero. 

Since L annihilates the vector with all entries 1, its 
determinant is zero. Kirchhoff's matrix-tree theorem Q 
and its generalizations express determinants of square 
submatrices of L as generating polynomials of spanning 
trees or rooted spanning forests in G. For any vertex 
i G V, let L(i) be the matrix obtained from L by deleting 
the ith row and column. Then Kirchhoff's theorem states 
that detL(i) is independent of i and equals 



det L(i) 
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where the sum runs over all spanning trees T in G. (We 
recall that a subgraph of G is called a tree if it is con- 
nected and contains no cycles, and is called spanning if its 
vertex set is exactly V.) The ^-independence of detL(i) 
expresses, in electrical-circuit language, that it is physi- 
cally irrelevant which vertex i is chosen to be "ground" . 
There are many different proofs of Kirchhoff's formula 
(JTJ ; one simple proof is based on the Cauchy-Binet the- 
orem in matrix theory (see e.g. 0)- 

More generally, for any sets of vertices /, J C V, let 
L(I\J) be the matrix obtained from L by deleting the 
columns / and the rows J; when I — J, we write simply 



2 



L{I). The "principal- minors matrix-tree theorem" reads 



det L(i±, . . . , i r ) = 



E 



n 



(2) 



where the sum runs over all spanning forests F in G com- 
posed of r disjoint trees, each of which contains exactly 
one of the "root" vertices i±,...,i r . This theorem can 
easily be derived by applying Kirchhoff 's theorem £[J to 
the graph in which the vertices ii, . . . , i r are contracted 
to a single vertex. Finally, the "all-minors matrix-tree 
theorem" (whose proof is more difficult, see ,2]) states 
that for any subsets /, J of the same cardinality r, 



detL(J|J) = e ( F ^ J ) II 

F£F(I\J) e£F 



(3) 



where the sum runs over all spanning forests F in G com- 
posed of r disjoint trees, each of which contains exactly 
one vertex from I and exactly one vertex (possibly the 
same one) from J; here e(F,I, J) = ±1 are signs whose 
precise definition is not needed here. 

Let us now introduce, at each vertex i € V, a pair 
of Grassmann variables ipi, r/>j. All of these variables 
are nilpotcnt (ipf — tpf = 0), anticommute, and obey 
the usual rules for Grassmann integration Writing 
T>(il), tp) = Yiiev dipi> we have, for any matrix A, 



T>{1>,$) e 



tliAt/> 



det A 



(4) 



and more generally 



■■ e(ix, . . .,v|ii, . . .,j r ) detA(n, . . .,i r \ji, ■ ■ ., j r ) (5) 



where the sign e(ii, . . . ,i r \ji, ■■•■>jr) = ±1 depends on 
how the vertices are ordered but is always +1 when 
(ii, . . . ,i r ) = (ji, ... ,j r ). These formulae allow us to 
rewrite the matrix-tree theorems in Grassmann form; for 
instance, (J2J becomes 



/ 



\a=l ) FeF(i!,...,i r ) e£F 

(6) 

Let us now introduce, for each connected (not necessar- 
ily spanning) subgraph T = (Vr,£"r) of G, the operator 



Qi 




n tiki** 



(7) 



family T = {Ti, . . . with I > 0, and let us try to 
evaluate an expression of the form 



It, e 



tpLtp 



(8) 



If the subgraphs Ti , . . . , T/ have one or more vertices 
in common, then this integral vanishes on account of 
the nilpotency of the Grassmann variables. If, by con- 
trast, the ri,...,rj are vertex-disjoint, then ex- 
presses / T>(t/j, f n'fc=l n ie Vr fc e ^ Llp as a sum 
over forests rooted at the vertices of Vr = (Jfe=i ^r* • 
In particular, all the edges of Er = Ui=i Er k must be 
absent from these forests, since otherwise two or more of 
the root vertices would lie in the same component (or one 
of the root vertices would be connected to itself by a loop 
edge). On the other hand, by adjoining the edges of Er, 
these forests can be put into one-to-one correspondence 
with what we shall call T-forests, namely, spanning sub- 
graphs H in G whose edge set contains Er and which, 
after deletion of the edges in Er , leaves a forest in which 
each tree component contains exactly one vertex from 
Vr- (Equivalently, a T- forest is a subgraph H with I con- 
nected components in which each component contains 
exactly one T^, and which does not contain any cycles 
other than those lying entirely within the Ti. Note, in 
particular, that a T-forest is a forest if and only if all the 
Ti are trees.) Furthermore, adjoining the edges of Er 
provides precisely the factor T\ 



eeE r 



u e . Therefore 



(9) 



where the sum runs over all T-forests H . 

We can now combine all the formulae © into a single 
generating function, by introducing a coupling constant 
tr for each connected subgraph T of G. Since 1+trQr = 



e *rQr we have 



2?(Vs$) e 



r = e (n*r) e 

r vertex- rer HeFr eeH 

disjoint 

(io) 

We can express this in another way by interchanging the 
summations over T and H . Consider an arbitrary span- 
ning subgraph H with connected components Hi , . . . , Hi; 
let us say that T marks Hi (denoted T -< Hi) in case Hi 
contains V and contains no cycles other than those lying 
entirely within T. Define the weight 



W(Hi 



E'r 



(11) 



(Note that each Qr is even and hence commutes with the 
entire Grassmann algebra.) Now consider an unordered 



Then saying that H is a T-forest is equivalent to saying 
that each of its components is marked by exactly one of 



3 



the r,; summing over the possible families Y, we obtain 
jV(1>,1>)e 

e (ri^Wjn^- ( i2 ) 

H spanning C G \»=1 / e£H 

H = (H u ...,Ht) 

This is our general combinatorial formula. Extensions al- 
lowing prefactors ipi 1 ipj 1 ■ ■ ■ tpi r tpj T are also easily derived. 

We shall discuss elsewhere some of the applications of 
(|12p. and restrict attention here to the special case in 
which tr — t whenever Y consists of a single vertex with 
no edges, tr — u whenever Y consists of two vertices 
linked by a single edge, and tr = otherwise. We have 

jT>(tp,^) exp tpLip + 1 ipiipi + u Wjj'ijjiipi'ipji/jj 

= e (nw^i+^i^i)) ( ia ) 

f £ f \i=l / eGF 

F= (Fi,...,F,) 

where the sum runs over spanning forests F in G 

with components Fi,...,F; here |VfJ and |FfJ 

are, respectively, the numbers of vertices and edges 
in the tree Fj. [We remark that the four-fermion 

term uJ2(ij) w ij4'i 1 Pi' l Pj' , Pj can equivalently be writ- 



where 

n ' = n + 5E 1 °g( 1 - T7r ") 

i 

= ^E'-^' • 7r -- - IE 77 ?" jE^'^K 2 

+ 0(^,7$. (14) 



When N = —1, the bosonic field 7r has —2 components, 
and so can be replaced by a fermion pair if we make 
the substitution 7Tj • 7Tj — > V'iV'j — iftiipj- Higher powers of 
Tvf vanish due to the nilpotence of the Grassmann fields, 
and we obtain the model if we identify t — — T, u = 
T. Note the reversed sign of the coupling: the spanning- 
forest model with positive weights (t > 0) corresponds to 
the cmtiferromagnetic iV- vector model (T < 0). 

An alternate mapping can be obtained by introducing 
at each site, in addition to the Grassmann fields tpi,4>i, 
an auxiliary one-component bosonic field (fi satisfying 
the constraint ipf + 2ttpiipi = 1. Solving this constraint 
yields ipi = 1 — tipiipi = e~ t ^ i ^ i and 

%J + 2t^j-l) = ^- Sfa - (1 - t^Vi)) 



(15) 



If we define the superfield ct, = (fi,tl)i,'^i) with inner 
ten, using nilpotency of the Grassmann variables, as product ct, ■ crj = tp%<pj + t{%jji%jjj — ipi'tpj), then the er- 



— (u/2)J2ij' l Pi'>PiLij'ipj'ipj-] If m = 0, this formula rep- 
resents vertex-weighted spanning forests as a massive 
fermionic free field 0, More interestingly, since 

\Vf, I — I-FfJ = 1 for each tree Fj, we can take u = —t 
and obtain the generating function of unrooted spanning 
forests with a weight t for each component. This is fur- 
thermore equivalent to giving each edge e a weight w e /t, 
and then multiplying by an overall prefactor fl^l This 
fermionic representation of unrooted spanning forests is 
the translation to generating functions and Grassmann 
variables of a little-known but important paper by Liu 
and Chow [ijj. 

The generating function of unrooted spanning forests 
is also of interest because it arises as a q — ► limit of the 
g-state Potts model, in which the couplings v e = e' 3Je — 1 
tend to zero with fixed ratios w e = v e /q [6j. 

Mapping onto Lattice a-Models. We now claim that 
the model (|13|1 with u — —t can be mapped, to all or- 
ders in perturbation theory, onto the ./V- vector model at 
N = — 1. Recall that the A 7 - vector model consists of spins 
lerJ = 1 located at the sites i £ V, with Boltz- 



mann weight e n where Ti = — T 1 ^ 



and T = temperature 
theory is obtained by writing cr 
with 7Ti G 



(ij) UJijfa-Uj 



1) 



pjV- 



Low-temperature perturbation 
1 and expanding in powers of n. Taking 



into account the Jacobian, the Boltzmann weight is e 



model with Hamiltonian ti = —T 1 X) (ij) w ij i&i ' &j ~ 1 ) 
and constraint er^ • ct, = 1 corresponds to the fermionic 
model (|13|l if we again make the identification t = —T, 
u = T . This CT-model, which is invariant under the su- 
pergroup 05P(1|2), has been studied previously by one 
of us [llj . It is presumably nonperturbatively equivalent 
to the N- vector model at N — —1, on the grounds that 
each fermion equals —1 boson. 

It is worth mentioning that the correspondence be- 
tween the spanning-forest model and these two er-models, 
while valid at all orders of perturbation theory, does not 
hold nonperturbatively. (This can be checked explicitly 
in the exact solution for the two-site model 12].) The 
error arises from neglecting the second square root when 
solving the constraints; we did not, in fact, parametrize 
a (super)sphere but rather a (super)hemisphere. Indeed, 
since t > corresponds to an antifcrromagnetic er-model, 
the terms we have neglected are actually dominant! But 
no matter: the perturbative correspondence is still cor- 
rect, and has the renormalization-group consequences 
discussed below. Furthermore, we conjecture that a non- 
perturbative correspondence can be obtained by using a 
tr- model with a suitable variant Boltzmann weight. 

Continuum Limit. Suppose now that the graph G is a 
regular two-dimensional lattice, with weight Wij = w > 
for each nearest-neighbor pair. We can then read off, 



4 



from known results on the iV- vector model |l3J, the RG 
flow for the spanning-forest model: it is 



dt 



2tt 



-P- 



(2tt): 



2.34278457- 



(2tt) 3 



-* 4 



1.43677- 



(2tt)< 



(16) 

where t — t/w and £ is the logarithm of the length rescal- 
ing factor; here the first two coefficients are universal (af- 
ter suitable normalization of the kinetic term), while the 
remaining coefficients are for the square lattice only. The 
positive coefficient of the P term indicates that for t > 
the model is perturbatively asymptotically free. Indeed, 
two-dimensional TV- vector models are asymptotically free 
for the usual sign of the coupling (T > 0) when N > 2, 
but for the reversed sign of the coupling (T < 0) when 
N < 2. Assuming that the asymptotic freedom holds 
also nonperturbatively, we conclude that for t > the 
model is attracted to the infinite-temperature fixed point 
at t = +oo, hence is massive and OSP( l|2)-symmetric. 
For t cr it < t < 0, by contrast, the model is attracted 
to the free-fermion fixed point at t = 0, and hence is 
massless with central charge c = —2, with the OSP(l\2) 
symmetry spontaneously broken. Finally, for t < i cr it we 
expect that the model will again be massive, with the 
05P(1|2) symmetry restored. 

More specifically, for t > it is predicted that the 
correlation length diverges for t | (or w f +oo) as 

1/3 



£ = a 



a (27r/3)(u>/i) 



1 - 0.0116221204- 



2tt w 

YT 

t 



0.00446142- 



(17) 



where is a nonperturbative constant (the terms in 
brackets are for the square lattice only). The numerical 
results of based on transfer matrices and finite-size 
scaling, are consistent with the nonperturbative validity 
of the asymptotic- freedom predictions l|16(l / i|17|l . but are 
inconclusive because the strip widths are small [L < 10). 
It would be interesting to make a Monte Carlo test of 
(|T7Jl . at large correlation lengths, along the lines of [H |. 

The numerics of are also consistent with the central 
charge c— —2 in the massless phase i crit < t < 0, but are 
not definitive because of the strong (1/log) corrections 
to scaling induced by the marginally irrelevant operator. 

Finally, the critical point t cr it presumably corresponds 
to the q — > limit of the antiferromagnetic critical 
curve in the g-state Potts model, under the identifica- 
tion w/t = [eP J — l)/q. Known exact results for the 
square lattice 0,0] yield (io/t) C rit = —1/4. The analy- 
sis of the critical theory proves rather difficult, but there 
are strong indications that it is simply a free OSP(\\2) 
model, i.e., the theory of a non-compact boson and a pair 
of fermions, with central charge c = — 1. 

Let us also remark that there exists a much-studied 
variant of the iV-vector model in which the high- 
temperature expansion on the lattice has been truncated 



so as to forbid loop crossings |17|. For — 2 < N < 2 
this model possesses several critical points; in particular, 
the dilute-loop critical point is expected to be generic in 
the sense that adding loop crossings acts as an irrelevant 
perturbation. For N = — 1 this yields a c — —3/5 theory 
|18j : the relation to the c = — 1 theory discussed above is 
mysterious and deserves further study. 

It would also be interesting to know whether our iden- 
tities are in any way related to the forest-root formula of 
Brydges and Imbrie 0], which leads to a dimensional- 
reduction formula for branched polymers. 
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